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We study two effective models for QCD, the Nambu-Jona-Lasinio -model and the linear sigma 
model extended by including a Polyakov loop potential, which is fitted to reproduce pure gauge 
theory thermodynamics, and a coupling between the chiral fields and the Polyakov loop. Thus 
the resulting models have as relevant degrees of freedom the Polyakov loop and chiral fields. By 
comparing the extended models with the bare chiral models we can conclude that the addition 
of the Polyakov loop is necessary in order to obtain both qualitatively and quantitatively correct 
results at finite temperatures. These results are extended to finite net quark densities, several 
thermodynamical quantites are investigated in detail and possible applications and consequences 
for relativistic heavy ion collision phenomenology are discussed. 



I. INTRODUCTION 

Phase transitions in strongly interacting matter have been a subject of intense theoretical, computational and 
experimental research over the past decades. Apart from the physical values of the QCD parameters, a lot of effort 
has been devoted to understand different limits of the theory as the mass parameters values, numbers of colors and 
flavors are varied. Several qualitative and in some cases quantitative aspects have been revealed: for example, we 
know that in the absence of quarks the SU(N) Yang-Mills theory has a global Zn symmetry [T|, and there exists a 
gauge invariant operator charged under Zn, the Polyakov loop, which can be identified as the order parameter of the 
theory. Hence, the deconfincmcnt phase transition can be easily characterized using the universality arguments [2]. 
Numerical studies have confirmed this picture as it has been found that the deconfinement phase transition is second 
order when the number of colors is iV c = 2 [3J, first order for -/V c = 3 [I] (although weakly [5]), and presumably first 
order for N c > 4 [B]. 

Adding quarks to the theory changes the picture considerably. For light fermions in the fundamental and pseudoreal 
representations for N c = 3 and iV c = 2, respectively, the corresponding Z 3 or Z 2 the center of the gauge group is not 
a good symmetry. However, the exact chiral symmetry of massless fermions is only little perturbed by small masses 
and the order parameter is the chiral condensate which characterizes the chiral phase transition. For iV c =3 and two 
massless quark flavors at finite temperature and zero net baryon density, the chiral phase transition is in the same 
universality class as the three dimensional 0(4) spin model [7], becoming a smooth crossover as small quark masses 
are accounted for [8]. For N c = 2 the relevant universality class is that of 0(6) [5]. 

In addition to finite temperature, one can also study the response of the QCD vacuum by considering finite net quark 
densities. At finite chemical potential one needs to take into account the pairing phenomena and superconducting 
phases result. However, a systematic approach based on the full QCD dynamics can be applied only in asymptotic 
densities |10] where asymptotic freedom simplifies the analysis. To determine the phases at intermediate densities 
relevant for phenomenological applications one needs to resort to effective models. Within the model studies taking 
into account the patterns of chiral symmetry breaking, one usually hopes to gain insight on the qualitative aspects 
of the phase diagram. Studies within different models have revealed that in cold dense matter, a first order phase 
transition to a superconducting phase characterized by nonzero diquark condensate takes place. Contrasting this 
with the finite temperature crossover transition at fx = one concludes that in the (T, /x)-plane there must exist 
a critical endpoint. A central paradigm for the first principle studies of the QCD equation of state [TT] over the 
(T, /x)-plane is therefore the existence, location and other properties of this critical point. The lattice determinations 
using different techniques have provided estimates for the location of the critical point. For two flavors see [TS] and 
[T3] . The existence of the critical point in three flavor QCD is currently debatable, [14] [15] , and in any case one should 
be careful in drawing any conclusions from Nf = 2 results to the physical 2+1 case . 

The perturbative calculations for cold quark matter at asymptotically large chemical potential or for hot quark gluon 
gas at high temperature cannot be directly applied to the phenomenologically relevant densities and temperatures. At 
finite temperature the perturbation expansion is known to converge poorly and reliable results can be obtained only 
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at near-zero chemical potential well above the suitably defined transition temperature T > 3T C , where the picture of 
weakly interacting dressed quark and gluon quasiparticles becomes correct. To address quantitative phenomenology for 
all temperatures at zero net quark density, a numerical " recipe" for interpolating smoothly between thermodynamics 
described by perturbation theory at high temperatures and by a resonance gas at low temperatures was proposed 
in |16| . taking also into account the contribution of quarks with finite masses. This approach provides one with a 
working definition of p(T, {m.;}), i = 1, . . . , Nf and of the resulting thermodynamics, but it does not yield insight to 
the nature of the underlying effective degrees of freedom and mechanisms responsible for the dynamics between the 
thermally dressed hot quarks and gluons and low temperature hadronic resonances. 

In this work we concentrate on the case of QCD with two light quark flavors and study a particular effective model 
description of it. This effective model is based on the following picture: As a function of the quark mass there exist 
two separate dynamical sectors. For massless quarks there is spontaneous breaking of exact chiral symmetry, while 
if the quarks are very heavy, and decouple from the dynamics, the center symmetry and its spontaneous breaking 
becomes relevant. In real QCD neither of these two is exact. Lattice calculations indicate that for quarks with 
finite masses, the transition is a smooth crossover as is also expected on the basis of universality arguments, and 
the transition can be located by measuring the expectation value of the chiral condensate. However, even if the 
discrete symmetry associated with deconfincmcnt is broken by the presence of light quarks, one can still study the 
temperature dependence of the Polyakov loop on the lattice. This has been done, and one observes the Polyakov loop 
to rise from zero to one as temperature is increased from low to high values. Due to this behavior one also speaks 
of deconfining phase transition [T7]. Moreover, the lattice results [T7] indicate that at zero chemical potential chiral 
symmetry breaking and confinement (i.e. a decrease of the Polyakov loop) occur at the same critical temperature. 
Several attempts to explain these behaviors exist 18J. Relying only on the exact and approximate symmetries of the 
system and general effective field theory methods a qualitative solution to this puzzle was established in [19] based 
on the idea of transfer of information from the order parameters to non-critical fields. As a function of quark mass, 
from light to heavy quarks, this mechanism allows one to treat either chiral symmetry or center symmetry as the 
relevant one driving the transition and through interactions allowing also the other would-be-order parameter field 
to behave in a similar way. The framework proposed in [19] also explains the independence of deconfinement and 
chiral symmetry restoration in the case of adjoint quarks which do not break the center symmetry. The (T, fi) phase 
diagram for adjoint two color QCD was considered in detain in [2"tJ] . 

This behavior has been reanalysed and confirmed in numerical studies of NJL model and linear sigma model coupled 
to the Polyakov loop via quarks which are integrated out in the random phase approximation [5TJ [22J [53] . In this 
work we study the coupling of Polyakov loop and chiral degrees of freedom described either with a linear sigma model 
or NJL model, and compare the resulting thermodynamics to that of two flavor QCD. Inclusion of the Polyakov 
loop makes the result qualitatively, and to some extent, even quantitatively insensitive to the underlying model used 
to describe the chiral degrees of freedom at finite temperature and small values of quark chemical potential. The 
importance of the Polyakov loop is not surprising, since the gluonic degrees of freedom are known to be important for 
the bulk thermodynamics of QCD matter when the net baryon densities are small. We consider also finite densities, 
and show that quantitative discrepancy between the two models increases as larger values of chemical potential are 
considered. To display this concretely, we determine the location of the critical point in (T, /x)-plane and show how 
the two models lead to very different results. 

The effective models studied here, may provide input to the phenomenology of relativistic heavy ion collisions. 
Recently, with the advent of RHIC data, it has been established that the spacetime evolution of the hot dense QCD 
matter is well described by nearly ideal hydrodynamics. This means in particular that the system evoles along the 
lines of constant S/N. Hence, we study in particular the behavior of isentropic lines in the (T, /i) phase diagram 
in these two models. Whether these have the tendency to focus on the critical point is important for the possible 
experimental discovery of the critical point in heavy ion collisions. We find that in these models strong focusing 
behavior does not exist. 

We introduce the models in some detail in section ITT] and carry out an analysis of the thermodynamics in section 



III Comparing these models at zero chemical potential, we find that they imply very similar results. We also compare 



with the results of the resummed perturbation theory [16] . At finite chemical potential we find that the quantitative 
results of the models show large deviations. Especially the location of the QCD critical point cannot be estimated 



reliably within these models. We end with concluding remarks and discussion of further prospects in section IV 



II. MODELS 



The chiral dynamics of two-flavor QCD is often formulated in terms of a linear or non-linear sigma model, which 
treats the Goldstone bosons as the relevant degrees of freedom. Of these two possibilities, the linear representation is 
more useful to study the finite temperatures and densities in order to find the phase diagrams of the theory, since also 
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the order parameter is included explicitly. Yet another possibility is to treat the fundamental fermion fields as basic 
degrees of freedom, the mesons appearing as the bound states of the theory, and this leads to Nambu Jona-Lasinio 
(NJL) -models. The effects of small quark masses are taken into account in the effective model Lagrangians by terms 
explicitly breaking chiral symmetry. These terms, appearing with small coefficients, render the chiral symmetry of 
the theory only approximate. Both of the above mentioned effective models for the phenomenology of two flavor QCD 
can be parametrized to describe equally well the vacuum structure at T = [i = 0. 

Finite temperature dynamics of SU (JV) pure gauge theory on the other hand is represented by Zjv symmetric 
effective theory for which the order parameter is the Polykov loop. Polyakov loop can be constructed and studied 
also in a theory with quarks eventhough the presence of fermions in the fundamental representation of the gauge 
group breaks the center symmetry explicitly due to the antiperiodic boundary conditions of the fermion fields at finite 
temperature. 

Therefore, in real QCD neither chiral symmetry or the center symmetry is exact, and we know that the finite 
temperature phase transition at /i = is a smooth crossover. However, one may ask which of the two symmetries is 
more accurate and would act as a "driving force" for the transition. Since chiral symmetry breaking is proportional 
to m qi and Zjq breaking is proportional to l/m q , in the case of two light flavors, it seems natural to consider the 
system to have an approximate chiral symmetry. This expectation is strenghtened also by looking at the spectrum 
of the bound states, as the pions clearly show the approximate Goldstone behavior. Based on these motivations, in 
[19j the situation was considered taking the m q = limit in which the chiral symmetry becomes exact, while the 
Zn symmetry is completely broken. Then the general principles of effective theory dictate a following form for the 
potential 

£[<7,7T<\$ = £ {<j,ir a } + £ o [0} + £ int [a,ir\0l (1) 

where £q[ct, 7r a ] is the chiral lagrangian which has exact chiral symmetry, £o[<p] is the potential for the Polyakov 
loop and contains both Zjq symmetric and symmetry violating terms and finally £ lnt is the part containing the 
interactions between the chiral fields and the Polyakov loop. As shown in [IS], the most important term for the 
dynamics is 4>(cr 2 + 7r 2 ), which leads to transfer of information between the order parameter and a non-order parameter 
field. When quark mass is increased away from the chiral limit, the transition becomes a smooth crossover, but the 
coincidence of the chiral symemtry restoration and deconfinement is expected as long as chiral symmetry remains 
good approximation. 

On the other hand one can consider infinitely heavy quarks, i.e. the pure gauge limit. Then the mechanism described 
above works similarly, but the roles of chiral symmetry and center symmetry are switched and the deconfinement 
order parameter drives the change of the chiral condensate. Again the two phenomena will coincide. Decreasing 
the quark mass from the pure gauge limit, the first order deconfinement line is expected to terminate at a critical 
point a some value of quark mass and for smaller values become a smooth crossover. Since lattice investigations find 
coincidence of critical temperatures related to chiral symmetry restoration and deconfinement for the accessible quark 
masses, it is reasonable to expect that there is a single phase border in (m q , T) plane interpolating between these well 
known small and large quark mass behaviors. The results of |19j can be applied to understand the behaviors near 
either small or large quark mass critical points. For intermediate values, more specific model studies or first principle 
lattice calculations are needed. 

In 21, 22 a specific model framework claimed applicable for all values of m q was proposed. In this work we study 
this framework in detail. We consider, side-by-side, both the NJL model and the linear sigma model (LSM) for two 
mass-degenerate quark flavors coupled to the Polyakov loop. The important feature underlying the dynamics in the 
approaches |21l 1221 123] is the assumption of independent deconfinement and chiral symmetry restoration described 
by the order parameters <p and a, respectively, and having independent effective potentials, U$ and U ai connected by 
interactions between the two. The central further assumption, then, is that the proposed interaction term yields the 
correct form for the resulting effective potential at all values of the quark mass, i.e. interpolates correctly between 
the limits of exact center symmetry and exact chiral symmetry hence describing also correctly the behaviors at the 
point (mphys,? 1 ), corresponding to real two-flavor QCD. We aim to test this underlying assumption in detail by 
cartographing the thermodynamics of these models over the (T, fi) plane and study in quantitative detail these two 
models against each other as well as against the numerical knowledge of real two-flavor QCD at zero net quark density. 
Let us start by describing the details of the models we use. To derive the grand canonical potential, we consider the 
following Lagrangian 



£ — -Cchiral + U<f,, (2) 

where we have separated the contributions of chiral degrees of freedom and the Polyakov loop. The part £ c hhai is 
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for the linear sigma model and NJL model, respectively, 

A 2 

£chirai = q{i^{d^^ig s A 8^)-g{G + i lb f-n))q~— {<r 2 + n 2 - v 2 ) 2 + H a , LSM (3) 
£chh-ai - q(*r(d^ - ig s A S^ ) - m Q )q - ^ NJL (4) 

where q = (u,d) is the light quark field, mo the bare quark mass (mo = m u = rrid, i.e exact isospinis assumed), a 
and 7? T = (it 1 , it 2 , it 3 ) constitute a chiral field E T = (a, if) and finally M = mo — G(qq). We work under the mean 
field approximation, hence the kinetic term of the chiral field is neglected in Q and the four-fermion interaction of 
the NJL-model lagrangian has been linearized in the condensate in Q. The symmetry breaking field H in ^ is 
H = /^m 2 , where /„■ = 0.093 GeV and m^ = 0.138 GeV. The coupling A 2 is determined by the tree level mass 
m 2 = 2A 2 / 2 + m 2 , which is set to be 0.60 GeV. In vacuum the expectation values of the fields are a = f„ and tt = 0. 
Requiring that the constitutent mass in vacuum is about 1/3 of the nucleon mass yields g — 3.3. In the NJL model 
the bare quark mass mo is taken to be 5.5 MeV and the coupling G — 10.08 GeV -2 . For a summary of the parameter 
values see table HI 

The Polyakov loop is included through the mean field potential 

U(cf>, <t>*,T)/T* = H 2 + <n + jM 2 ) 2 , (5) 



where 



6 2 (r)^ ao + «iY+«2(^) 2 + «3(^) 3 , (6) 



and the constants cii,bi are fixed to reproduce pure gauge theory thermodynamics with phase transition at To = 270 
MeV. We adopt the values determined in [25] and shown for completeness in table [i] Here <ft is the gauge invariant 
Polyakov loop in the fundamental representation. One could also include other loop degrees of freedom, say, adjoint 
or the sextet. Here we choose to start with the mean field potential of the fundamental loop parametrized to describe 
the pure gauge thermodynamics and study how the interactions with the chiral degrees of freedom affect it and 
compare against full QCD thermodynamics. The extensions towards other possible degrees of freedom we leave 
for future work. Recently there has been interesting developemcnts in the construction of effective theory for the 
pure gauge thermodynamics, e.g. [24j[25], which could be in principle coupled to chiral fields to obtain an effective 
theory for QCD. Here we choose to remain with the mean field potential §5§ which seems to capture the pure gauge 
thermodynamics sufficiently well for our purposes. 

Then, for a spatially uniform system in thermodynamical equilibrium at temperature T and quark chemical potential 
fi the partition function is 

Z = Trcxp[-(W- M A0] 

= JvqDqexp[J(^ + m°q)]- (7) 

The integration over the spacetime in the action is over the compact euclidean time direction and over the spatial 
three- volume V. Since the action is quadratic in quark fields, the functional integral is easily performed with standard 
methods leading to the grand canonical potential 



T\uZ 

t^chiral + U ^ + Qq q , (8) 



where 



[/chiral = ^(c7 2 + n 2 -v 2 ) 2 -Ha, LSM 

C/c hiral = { ^^, NJL (9) 

for the chiral contribution and U<p = U(<fr, <fi*,T). The interactions between the Polyakov loop and chiral degrees of 
freedom are 



n qq = -2N f T J ^( T r c ln[l + L e -(^)/ T ]+Tr c ln[l + Lt e -(^)/T^ _ QNf j ^E8(A 



|pf ), (io) 
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0.093 GeV 


4.44 


3.3 
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5.5 MeV 


651 MeV 


10.08 (Gev)" 2 
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ai 


0.1 




6.75 


-1.95 


2.625 
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-7.44 


0.75 


7.5 



TABLE I: The parameters used for the effective potential 



where the trace over color remains, E = y/p*+M 2 (and further M — go in LSM). In LSM we neglect the vacuum 
contribution term in £lq q and in the NJL model it is controlled by the cutoff A as indicated in (10 1. Performing the 
remaining trace gives 

(Tr c ln[l + Le-( E -^' T ] + Tr c ln[l + Lt e -( B +^]) = 

ln(l + 3(0 + 0* e -(E-*O/T) e -(E-,O/T + e -3(E-^/ T) + m(1 + + ^-(E+^/Ty-iE+^/T + ^{E+^/Ty (n) 

Note that in principle the chemical potential affects the Polyakov loop potential directly, see [23], but we will not 
be considing these effects. Having determined the grand canonical potential both for the Polyakov loop linear sigma 
model (PLSM) and Polykov loop NJL model (PNJL) , the thermodynamics is now determined by solving the equations 
of motion for the mean fields, 

on an on 

^= ' ^r°< ^= ' < 12 > 

and then the pressure is given by evaluating the potential on the minimum: p = —Q(T, /i). We now proceed to solve 
numerically the thermodynamics of PLSM and PNJL models and compare them with each other as well as against 
the numerical results on two-flavor QCD at zero chemical potential. 



III. NUMERICAL RESULTS 



A. Thermodynamics at /i = 0, comparison to QCD 

Let us first consider the models at zero chemical potential but finite temperature. When finite temperature is 
considered, it is well known that the quantitive results of chiral effective theories differ: While NJL model predicts 
chiral restoration at T ~ 150 MeV, the linear sigma model leads to result T ~ 190 MeV [8]; allowing for finite 
chemical potentials only widens the spread. Including the Polyakov loop has important consequence as now both 
of these models predict a crossover near T ~ 210 . . . 230 MeV within 20 MeV of each other, as we show in Fig. [TJ 
where we plot the temperature derivatives of the condensates. The location of the peak defines a critical temperature 
around which the crossover takes place. 

Then consider the pressure. The coupled models, PNJL and PLSM, have already been shown to agree with lattice 
data at and above T c fairly well [22J [23] with some fine-tuning. Namely, lattice data implies a critial temperature 
T c ~ 175 MeV, a value ~ 20% lower than we obtain. In [22] it has been noted that better agreement can be achieved 
in PNJL model by detuning the Polyakov loop potential away from the pure gauge thermodynamics through shifting 
of the parameter T down to 190 MeV. Since the lattice data has still some uncertainty to it due to the extrapolation 
to continuum limit, we choose not to aim for perfect fits and rather plot the thermodynamical quantities as a function 
of T/T c when comparing different models. Actually, the value of T c should be determined by allowing for additional 
degrees of freedom below T c not considered in these effective models. We have not included the finite temperature 
contributions of the pions and more massive resonances, since we do not have a dynamical way to decouple their 
contribution at high temperatures. However, our aim is to study the interplay of chiral fields and Polyakov loop 
and the resulting thermodynamics for the temperature range T c < T < 3T C for which the coupled models under 
consideration seem to work well. 

Here, to compare with QCD, we use the numerical result of |16j in which the high temperature part of the curve is 
based on the full 0{g 2 log g) calculation in pure gauge theory, supplemented with a more phenomenological "recipe" 
to include the contribution of Nf massive quarks at order 0(g 2 ). At low temperatures this result is matched smoothly 
on the resonance gas result from [25]. This QCD+resonance gas result is shown in the left panel of Fig. [2] at Nf = 2 
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FIG. 1: Left Panel: Temperature derivatives of the mean fields 4>{T) and cr(T) at fi — in the PLSM model. Right Panel: 
Same observables as in the left panel but in the PNJL model. 



appropriate for this work |34j . Also shown are our results for the pressure obtained from the coupled models PLSM 
and PNJL as a function of t = T/T c . For comparison we show the corresponding results from LSM and NJL models 
without the Polyakov loop. Observe how the inclusion of the Polyakov loop increases the result for the absolute value 
of the pressure by roughly 80 %. The addition of the Polyakov loop is therefore necessary in order to quantitatively 
obtain the required rise in the pressure towards the Stcfan-Boltzmann limit of QCD corresponding to the horizontal 
dashed line in the figure. The difference between the chiral model pressure and the additional increase due to the 
Polyakov loop can be understood by looking at the contribution of bosons and fermions in the ideal gas result which 
for zero chemical potentials is 

71-2714 f,„2 ^ 7N c N f \ 

Setting N c = 3 and Nf = 2, the ratio of the bosonic and fermionic contributions is QbI 9f — 16/21 w 0.76. 

Both PNJL and PLSM models give a good overall desctiption of QCD pressure above T c . Let us then turn to the 
analysis of more differential observables. The entropy density s(T) = p'(T), the energy density e(T) = Ts(T) — p(T), 
trace of the energy momentum tensor O^JT) — e(T) — 3p(T) and the heat capacity c(T) = e'(T) = Tp"(T), to name 
a few, serve as important probes which, together with lattice data, allow to probe the validity of different models and 
the knowledge of the temperature dependence of these quantites is important for phenomenology of rclativistic heavy 
ion collisions as well as for cosmology. 

Let us first consider the trace of the energy momentum tensor. We show the result of different models in right panel 
of Fig. [2] Notice how also here the inclusion of the Polyakov loop is imperative to match lattice data in comparison 
to LSM and NJL models: First, the addition of the Polyakov loop removes the qualitatively different structures in 
LSM and NJL model around T c and in the coupled models a smooth universal curve results. Second, the inclusion 
of the Polyakov loop is important also for obtaining agreement with the asymptotic behavior above T c . In fact both 
LSM and NJL models alone predict that at high temperatures O^JT) x T 4 = const, contrary to the observation. 

The information contained in the trace anomaly can also be represented in terms of the " equation-of-state - 
parameter" w(T) = p(T)/e(T) which we show in Fig. [3] Both PLSM and PNJL models give again very similar 
results, significantly lower than the pQCD result also shown in the figure. The effective model results are similar to 
recent lattice data [27 . There are two clear features deserving further numerical and theoretical studies: First, the 
drop near T c leads to w ~ 0.1 in PNJL and LSM models as well as in the lattice data while the perturbation theory 
result is larger, w ~ 0.15. Second, below T c the pQCD model leads to larger rise in w(T) than the lattice data and 
PNJL and PLSM models. These features arc likely to be affected by the resonance gas dynamics neglected in the 
effective models but which are present in the pQCD model. 

The information contained in the derivatives defined above can be presented and explored in various ways. We 
choose to follow the presentation in [16 , since that allows also a quantitative comparison against the results from the 
resummed perturbation theory. In Fig. [4] we plot the effective numbers of degrees of freedom defined by 

„ - < T ) h - S ( T ) ,• - < T ) ™ 

9eS — [ir 2 T 4 ~\ ' cff — [27r 2 T 3 ] ' cf ^ — |"2jr 2 T 3 l l 1 ^/ 
L 30 J L 45 J L 15 J 
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FIG. 2: Left panel: Pressure from the models with and witout the Polyakov loop at fi — 0. Also shown is the curve interpolating 
between the resonance gas and resummed perturbation theory results as well as the constant corresponding to the Stefan- 
Boltzmann limit of two-flavor QCD. Right panel: Similar figure for the trace anomaly (e — 3P)/T 4 . 
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FIG. 3: The equation of state parameter w(T) =p(T)/e{T). 

Looking at the effective degrees of freedom, we again see that both PLSM and PNJL model results are consistent 
with each other and with the corresponding QCD+resonance gas results. Only i c g shows some qualitative differences 
in how these two models respond to finite temperature as the peak at T c is sharper in PLSM model than in PNJL 
model. This peak arises since i e g is proportional to heat capacity which diverges in a second order phase transition. 
Another issue present in the temperature dependence of i e g is a small second peak visible very weakly in PNJL 
model but more strongly in PLSM model. This is due to the remnant of the deconfinement transition described 
by the temperature dependence of the parameters in the Polyakov loop potential. Note how it is vital to look at 
more differential observables, second derivatives in this case, to see this effect. More precise lattice data is needed to 
determine if the structure of two independent phase transitions connected by interactions underlying these effective 
models is indeed correct. 



B. Thermodynamics at /i 7^ 0: location of the critical point? 

Let us then consider the consequences of nonzero net quark density by allowing for finite chemical potential. As a 
starting point we will use the grand canonical potential derived in section [TTJ and neglect possible direct ^-dependence 
of the Polyakov loop potential which has been discussed eg. in ref. [23 . 
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FIG. 4: Left Panel: g e «{T) at /x = 0. Middle Panel: h cS (T) at fi = 0. Right Panel: i eS (T) at /j, = 0. 



First we study if the coincidence of deconfinement and chiral symmetry restoration holds also at finite chemical 
potential. In Fig. [5] we show the derivatives of the condensates as a function of temperature at ^ = 100 MeV and 
in Fig. [6] at /i = 250 MeV. We observe a coincidence of peaks in the derivatives of the condensates in both of these 
effective models. Note that at finite chemical potential <f> and (jf are no more equal. At large chemical potentials 
we observe that already at first derivatives of the condensates a double peak structure arises. This is due to the 
fact that as chemical potential is increased, the critical temperature in the chiral sector decreases as shown by the 
location of the leftmost peak in the figures, while the remnant of the deconfinement transition in the Polyakov loop 
potential is unaffected by the value of the chemical potential and remains visible at temperature T ~ 200 MeV. As the 
two distinct transitions underlying these effective models are separated over wider temperature range their separate 
features become also more visible. Hence, this provides another way to numerically investigate the correctness of the 
initial assumption of independent deconfinement and chiral symmetry restoration underlying these models. 




FIG. 5: Left Panel: Temperature derivatives of the mean fields 4>{T) and a(T) at \i = 100 MeV in the PLSM model. Right 
Panel: Same observables as in the left panel but in the PNJL model. 

We show the (T, \i) phase diagram in Fig. [7j where the solid line shows the result of PLSM and the dotted line 
shows the result of PNJL model. The lines have been determined by finding the location of the peak in the chiral 
condensate a(T) at each chemical potential. Low temperature and low net quark density phase is confined and chiral 
symmetry is broken, while at high temperature the chiral symmetry is restored and the system is deconfined. The 
two models yield very different values for the location of the critical point: while the PLSM gives (T c , /x c ) = (195, 141) 
for the position of the critical point, the PNJL model yields (88,329). This can be due to several reasons, e.g. the 
neglect of the possible chemical potential dependence of the Polyakov loop potential. On the other hand, this can be 
due to neglect of possible relevant degrees of freedom. It is likely that at finite chemical potential diquark degrees 
of freedom become important and should be taken into account; see |28j . This situation would be similar to the 
differences between LSM and NJL models at finite temperature which were shown to reduce once the Polyakov loop 
dynamics is accounted for. 

To discover the critical point experimentally in heavy ion collisions, it would be desirable to have a reliable quanti- 
tative theoretical estimate of its location. As we have explicitly seen here, various effective theory estimates deviate 
a lot when finite chemical potentials are considered. The lattice determinations using different techniques also lead 
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FIG. 6: Left Panel: Temperature derivatives of the mean fields (f>(T) and a(T) at fi — 250 MeV in the PLSM model. Right 
Panel: Same observables as in the left panel but in the PNJL model. 
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FIG. 7: The (T, n) phase diagram. Leftmost pair of curves shows the result for LSM and PLSM while the second pair is for 
NJL and PNJL models. The solid part of the curve denotes a first order transition while the dashed part is a crossover. 



to very different results for the location of the critical point. For example, for two flavors the authors of [T2] find the 
critical point at hb ~ 360 MeV while in [13] a value \ib ~ 180 MeV is reported. Currently the existence of the critical 
point is debatable |14l 115) . and in any case one should be careful in drawing any conclusions from Nf = 2 results to 
the physical 2+1 case . 

With these remarks in mind, let us assume that the critical point in the (T, f£) plane exists as implied by these 
effective theories. Then, even if the exact location of the critical point is not exactly known, one may argue in favor of 
its experimental detection if the spacetime evolution of the strongly interacting elementary particle matter is such that 
the system passes through the vicinity of the critical point starting from almost any initial condition. The outcomes 
of such focusing behavior have been recently advocated for in [29J strongly motivated by [30 . However, the focusing 
observed in |30j can be due to the particular equation of state applied in that work and is not a general feature of 
hydrodynamics approach applied successfully to describe the RHIC data as in e.g. [31]. Given the success of ideal 
fluid hydrodynamics in the description of the RHIC data, it is likely that the system expands nearly isentropically. 
Therefore, to decide if the focusing behavior in the models studied here should occur, we find the adiabats of PLSM 
and PNJL models in the (T,/j) plane. The result is shown in Figs 8] and [9] and the conclusion for both of these 
models is negative. Based on these figures we conclude that for the hydrodynamical evolution starting from a near 
zero net quark density, as would be the case in a central Au+Au collision at RHIC, the trajectory of the system in 
(T, /z)-plane is non-focusing. As the closeups in the right panels of Figs [8] and [9] show, there is no special behavior 
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near the critical points. This result is similar to the one obtained on the lattice [27] , 
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FIG. 8: Left Panel: Constant S/A curves in PLSM model. Right Panel: A closeup on the critical point at (T c ,/i c ) = (195, 141). 
Every second curve has been drawn with the dashed line only to enhance readability. 
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FIG. 9: Left Panel: Constant S/A curves in PNJL model. Right Panel: A closeup on the critical point at (T c ,fj, c ) — (88,329). 

Another interesting quantity is the sound speed shown in Fig. [l0]for chemical potentials near the critical one. The 
stronger dependence of the location of the minimum in c s (T) as a function of the chemical potential in PNJL model 
is due to the fact that near the critical point the value of T c (/i) changes more rapidly than in the case of PLSM 
model. Together with Figs. [8] and [9] we see from here that in the hydrodynamical evolution the sound speed is very 
small during the part of the evolution the system spends near the phase transition region. At (T, /i) regions relevant 
for RHIC and LHC/ ALICE phenomenology, the sound speed is probably well approximated by the [i — result as 
the system evolves very close to the n ~ axis along a nearly parallel trajectory until low temperatures deep in the 
hadronic phase are reached and the trajectory bends to end at the finite value /i vac at zero temperature. However, if 
it was possible to create systems which would follow the trajectory bending along the phase boundary as the ones at 
large chemical potential in Figs [8] and [9j then the sound speed could stay small over larger temperature range and this 
might have consequences for e.g. Mach cones created by high momentum jets traversing the thermal medium [32J. 
On the other hand, while the collisions planned at GSI/FAIR facility might lead to more optimal S/N trajectories 
for this argument to work, the probability for the production of the required high momentum probes is much smaller. 
However, we remind the reader that a thorough investigation of the heavy ion phenomenology should be carried out 
with an effective theory suited to describe the Nf = 2 + 1 case. 
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FIG. 10: Left Panel: Sound speed at different chemical potentials in PLSM model. Right Panel: Same for the PNJL model. 



IV. CONCLUSIONS AND OUTLOOK 



In this paper we have analysed thermodynamics of two effective theories for QCD containing as relevant degrees of 
freedom the Polyakov loop and chiral fields within a framework proposed to correctly interpolate between the pure 
gauge, center symmetric, and chirally symmetric two flavor QCD. At zero chemical potential but finite temperature 
we observed that the contribution from the Polyakov loop to thermodynamical quantities like the pressure is very 
important; roughly half of the total pressure. Another aspect which underlines the importance of the Polyakov loop 
sector is the fact that adding it to the dynamics tends to diminish the qualitative differences present in the bare 
chiral theories at finite temperature. We have considered explicitly two realizations, the linear sigma model and the 
NJL model. Our results indiate that even if both LSM and NJL models describe the chiral dynamics of the QCD 
vacuum correctly, they must be supplied by other degrees of freedom in order to obtain quantitatively correct effective 
description of QCD thermodynamics. We presented a comparison with respect to a recent result interpolating between 
resummed perturbation theory result and resonance gas result |16j , and outlined how future lattice simulations could 
allow one to obtain more insight into the QCD dynamics near the phase transition. 

We also considered finite chemical potentials and determined the (/x, T)-phase diagram of these theories. Here 
we observed large discrepancies between the PLSM and PNJL models. Based on the above discussion, an obvious 
explanation would be that some important dynamics is again being missed. A natural candidate for a new relevant 
degree of freedom is the diquark condensate responsible for the color superconducting phenomena. We aim to extend 
our work towards this direction next. We also determined the lines of constant S/A in (T, /i)-plane and discussed 
the implications for the equation of state as well as for the possible focusing behavior relevant for the experimental 
discovery of the QCD critical point. 

These results can be used for phenomenological applications. The parametrization of the equation of state obtained 
from these models has been shown to agree with lattice data at /i = but is easily evaluated also at finite values of the 
chemical potential. Hence it could be applied in hydrodynamical simulations of ultrarclativistic heavy ion collisions. 
We have evaluated the sound speed, and shown that on the hydrodynamically relevant trajectories there may be 
substantial temperature range over which the sound speed is small. However, it may prove difficult in the laboratory 
to create systems which would follow these particular trajectories. We hope to extend these phenomenological studies 
within realistic hydrodynamics in near future. 

There are several improvements to be addressed. We have already stressed that below the critical temperature, a 
more careful treatment of hadronic degrees of freedom is required, and the lack of this treatment is best seen in the 
apparent underestimate of the pressure below T c . Another issue clearly concerns the number of active flavors. Here 
we have concentrated only on the case of two flavors, while for more quantitative phenomenology it is vital to have the 
effects of the strange quark under control. However, as a first approximation for the equation of state at small chemical 
potential it should be reasonable to multiply the pressure by the overall factor gsB(Nf = 3)/gsB(Nf = 2) where g 
counts both the bosonic and fermionc degrees of freedom. Lattice data [33] indicates that the flavor dependence of 
the QCD pressure is dominated by the Stefan-Boltzman factor and the above scaling can be applied within 10 % 
accuracy. Similar reasoning has been used also in [TB] to obtain the perturbative QCD result for full Nf flavors with 
masses {tojv> }. We aim to address these issues in future work within the effective theory framework discussed here. 
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